Selection of the Taylor-Saffman Bubble without Surface Tension 
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By presenting the most general to date class of exact solutions, we demonstrate that the attractor 
for zero surface tension non-singular bubble dynamics in a Hele-Shaw cell is the selected pattern. 



Introduction. Steady patterns, often observed in var- 
ious unstable non-equilibrium systems, are formed, as a 
rule, after a transient nonlinear interface dynamics be- 
tween two phases. The problem of selecting a single ob- 
served pattern from a continuum of possible solutions, 
such as a parabolic dendrite in supercooled liquid [l[ and 
a finger in a Hclc-Shaw cell @, attracted considerable at- 
tention in the past (see details in Q ) . Besides the obvious 
importance for non-equilibrium physics, the interest was 
heated up by mathematical difficulties to include surface 
tension, considered indispensable for the selection. These 
difficulties were resolved in Q , and it was shown that sur- 
face tension indeed selects the pattern. More recently, 
analysis of exact solutions for the interface dynamics in 
a Hclc-Shaw channel demonstrated that selection of the 
observed finger Q is possible without surface tension [j| . 
The present work significantly extends this result by ad- 
dressing the bubble dynamics instead of a finger. 

The pattern selection problem for an inviscid bubble 
dragged by a viscous flow in a narrow gap between par- 
allel plates (Hele-Shaw cell) was posed by Taylor and 
Saffman in 1959 @ and was later addressed experimen- 
tally [3] and theoretically To be specific, the problem 
was how to select, from the continuum of steady state 
solutions obtained for zero surface tension, the unique 
bubble Q with velocity, twice the background flow veloc- 
ity in the limit of low surface tension. While it has long 
been known that the inclusion of surface tension leads to 
velocity selection Q, we demonstrate in this paper that 
the selection mechanism does not require surface tension. 
Instead, the selection comes about because the selected 
pattern is the only stable fixed point (attractor) of the 
non-singular bubble dynamics without surface tension. 

Problem formulation and plan of the paper. A top- view 
of a Hele-Shaw channel with lateral walls at y = ±7r in 
our scaled units and with the bubble moving to the right 
is shown on Fig. la. The fluid (oil) velocity obeys the 
2D Darcy law, v = — Vp, where p is scaled pressure. 
Far from the bubble the oil flows along the Ox axis with 
uniform velocity, V = 1, thus p = — x, when \x\ — > oo. 
Owing to incompressibility, V • v = 0, p is harmonic, 
V 2 p = 0, in the viscous domain, D(t), where t denotes 
time. It is thus convenient to introduce a complex po- 
tential, W = $ + i^, where <f> = — p and \I/ is the stream 
function. In view of the uniform far-field velocity, one has 



W = z for \x\ — > oo, so ^ = ±7r at y = ±ir, since the lat- 
eral walls are streamlines. Because pressure is constant 
(taken to be zero) in the inviscid bubble and continuous 
across the oil/bubble interface, T(t), if surface tension 
is neglected, then $ = at T(t). The fluid domain in 
the VF-plane is shown in Fig. lb, where the vertical slit 
maps to the interface, T(i), and the rest of the horizontal 
strip, — 7r < ^ < 7r, to the domain, D(t), in the physical 
plane (Fig. la). The kinematic identity, V n = v n , stating 
the equality of the normal velocities of the interface, V n , 
and of the fluid, v n , completes the formulation of this 
free-boundary problem of finding T(t) given T(0). 

This long-standing nonlinear unstable problem [9| was 
shown to have an integrable structure and to possess 
deep connections with other branches of mathematical 
physics (Toj . While numerous exact solutions were ob- 
tained for time-dependent interfaces (listed in [ll[) and 
for steady shapes (see [l2| and references therein) , almost 
all of them were for singly-connected domains D(t). 

Here we present exact solutions for an evolving bubble 
in a Hele-Shaw cell when the fluid domain D(t) is doubly 
connected (see Fig. la). We then show that these solu- 
tions explain not only how the moving bubble eventually 
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FIG. 1. The fluid domain D(t) for a moving bubble in a 
Hele-Shaw channel (a), and the corresponding domains in the 
complex potential plane (b) and in the auxiliary £-plane (c). 
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reaches a stationary shape when t — >• oo, but also w/iy 
the selected bubble moves precisely twice as fast (if to 
neglect surface tension) as the background flow Q ■ 

In what follows, we first address the dynamics of a bub- 
ble with a shape symmetric with respect to the channel 
centerline. The exact solutions we obtain for this sim- 
pler case illustrate our main result, namely, that all non- 
singular solutions describing bubble shapes converge to 
the only stable fixed point of this infinitely-dimensional 
dynamical system, which is precisely the selected mem- 
ber of the continuous family Q. After that, we extend 
our solutions to asymmetric bubbles and arrive at the 
conclusion that U — 2 is still the selected value. 

The conformal map and the equation of motion. For 
a symmetric bubble we introduce a conformal map z = 
f(t, C) from the exterior of the unit circle with a cut in 
the complex plane, £ = £ + (see Fig. lc), onto the fluid 
domain, D(t) in the z-plane, z = x + iy. The unit circle, 
|C| = 1, maps onto the interface, T(t), and the cut sides, 
C = £ ± iO, where 1 < l/v(t) < £ < oo, map onto the 
channel walls, y = ±7r, so that C = 1/V and £ = +oo are 
mapped to x = +oo and x = — oo, respectively. Thus, the 
polar angle (ft = arg(£) parametrizes T(t): z = f(t,e % ^). 
It is easy to see that the complex potential, W = $ + i^, 
satisfying all aforementioned boundary conditions, is 

^(i, C )=logl^. (1) 

Rewriting the normal velocity of the interface, V n , as 

V„ = V x h - Vih = Xm(Vl) = Xm(z t z 3 ), 

where I = l\ + U2 = dz/ds is the unit tangent vector 
along the interface, T(t) = z(t, s), parametrized by an 
arclength, s, and the subscripts stand for partial deriva- 
tives. The normal fluid velocity we rewrite as 

v n = <9„$ = <9 S *, 

by virtue of the Cauchy-Riemann condition. Equating 
the two last formulae, as required by the kinematic iden- 
tity, and reparametrizing s — > (ft, we obtain Xm{ztZ^) = 
Calculating ^ from (1) for £ = e 4 ^, we obtain the 
equation for the moving interface, z(t, (ft) — f(t, e 1 ^): 

Xm{z t Zs) = He— — . (2) 

For stationary solutions, z(t, (ft) = Ut + Z((ft), where the 
velocity, U, is a constant, equation (2) is simplified: 

Xm(^) = ^Tle-^—. (3) 

The solution of (3) is the sum of two logarithms: 

Z = -log(l - ve*+) + obg(l - ve-*+), (4) 

where the coefficient at the first term equals 1 to satisfy 
W = z in (1) when z — > 00. Expression (4) is precisely 



the one-parametric continuum of stationary bubbles ob- 
tained in Substituting (4) into (3), we obtain 

1 + a 

We will show that all solutions with a/0 are unstable, 
and, if perturbed, move to the solution with a = 0. This 
will solve the selection of U = 2 in accordance with (5). 

The finite-parametric solution. Being integrable, equa- 
tion (2) possesses a rich list of exact solutions, and we 
found them. Many of these solutions blow up in finite 
time. Leaving those aside as physically non-realizable, we 
present here a rich class of finite-parametric non-singular 
solutions (analogous of those obtained earlier [l3[ for fin- 
ger dynamics), which remain finite for all times: 

jv 

z = r(t) - log(l - i/(t)e**) + «fe lo g(! - a k {t)e~ 1 *), 

k=0 

(6) 

where ao — ct and r are both real, cto = v, \ak\ < 1 for all 
times, and otk are constants. These parameters must be 
chosen so that critical points of the conformal map always 
stay inside the unit circle to prevent blow ups [lj]. Also, 
the symmetry of the bubble requires that for each term 
with complex ctk and aj, in the sum there is another term 
with otk and et^. It is easy to verify that (6) is indeed a 
solution of (2), where the time-dependence of r, v, and 
Ok is given by the following set of N + 2 equations: 

N 

Pk = r + log d k /{a k - v) +^a ; log(l - a ; a fe ), (7a) 

1=0 

N 

2t + 2t = (1 + a)T - log(l - i/ 2 ) + J2 "fc lo g(! - ™ k ), 

k=l 

(7b) 

N N 

A = (H 2 - 1) log(l - v 2 ) + ak ® 1 ~ afeS ')' 

fe=l 1=1 

(7c) 

where k = 1, ...,7V in (|7a|) . Here the f3k, the initial time, 
to, and the bubble area, A, are the constants of motion. 

The attractor. It follows from (7b) and (7c) that r(t) 
diverges linearly with time, t — > Ut, when t — > +00. 
Since /3k is a constant, a real part of at least one logarithm 
in (7a) should go to —00 in large times so as to compen- 
sate a divergent positive r(t). This is possible only if all 
afc(t) — > for k > 0, as t — > +00. Thus, we conclude that 
the origin, £ = attracts all auit) different from v(t). 
But in this case the solution (6) approaches the family of 
stationary bubbles (4) discussed above, namely 

Of 

z = logfl - ve 1 ^) +cdog(l - ve~^), 

I + a 

which describes the selected value, U = 2, only if a = 0. 

If ao(0) = i/(0) then ao(t) = v{t) for all times. To test 
the stability of the trajectory, ao(t) = v(t), one has to 
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deviate a little ao from v. One can readily check that the 
dynamics of ao, v and r is given in this case by 

Po = T + loga /(a -i/)+cdog[(l- |a | 2 )(l - a 2 ,)], 
2(t + t ) = r - log(l - j/ 2 ) + a log |1 - a i/|, 
A = - log(l - ^ 2 ) + (a 2 /2) log[(l - |a | 2 )|l - ag|]. 

These equations clearly show that ao — > 0, r — >■ 2i, and 
i/ — > \/l — when t — > oo, so the interface becomes 



z = 2t - log(l - i/e^). 



(8) 



This is precisely the selected pattern with U = 2. Thus, 
£ = v repels nearly located singularities, which move to- 
ward 0. Therefore, the selected bubble, described by (8), 
represents the only attractor, £ = 0, of the non-singular 
subset of the finite-dimensional dynamical systems (6). 

The infinitely-parametric solution. Let us extend (6) 
to the case of infinite number of parameters: 



r(t)+ log 



1 



1 - v{t)e 



A f") p k (t,w) 



fc=0 



where ao = v and \a^\ < 1 for all k. Solution (9) coincides 
with (6), if the functions pk(t,w) are constants, and con- 
tains all previously known solutions for single-connected 
and symmetric double-connected cases. (Apparently one 
can approximate any non-singular solution of (2) by (9), 
although we have not proved this.) The constants of mo- 
tion for (9), analogous to p k in (7), are B k = f(l/a k (t)), 
which yield 



B k = r(t)+log 



Ok(t) 



1 -!/(*)**(*) f ^ 



ai(t) Pi{t,w) 



w - l/ojfe(i) 



for A; = 1, N. A proof, that r diverges linearly with t, 
as t — > oo, and testing the points, v and 0, for stability is 
the same as above with the same conclusion that all a& — > 
0, implying that an arbitrary shape, expressed by (9), 
moves toward the selected bubble (8) with U = 2. After 
having analyzed symmetric bubbles, we now move on to 
the non-singular dynamics of non-symmetric bubbles. 

Non-symmetric case. The complex potential for a non- 
symmetric shape requires infinity of reflected images, so 
we map conformally the annulus, < Jq < \(\ < 1, in 
the £-plane (see Fig. [2j) onto the fluid domain, D(t), in the 
z-plane, so that the inner circle, |£| = y /q 1 is mapped onto 
the interface T(t). The unit circle, \Q = 1, is mapped 
onto the channel walls, y = tt and y = 0. Under reflec- 
tion with respect to the unit circle (see Fig. 2), we obtain 
the pre- image of the domain, D(t), which is the complex 
conjugate of the D(t). We map 

{Vq < ICI < i/Ve} -> {D(t) u D(t)}, 

where the annulus, ^/q < |£| < 1/y/q is cut along the 
part of the unit circle, < arg £ < 7, so that the inner 
(outer) cut side is mapped onto the north wall (its mirror 




FIG. 2. Flow domain (shaded region) in the £-plane for an 
asymmetric bubble; see text. 



image), where y = \& = ±7r, while the complimentary arc 
along the unit circle, 7 < arg£ < 2tt, is mapped onto the 
south wall, where y = 'J = 0. We fix the map by sending 
the points £ = 1 to x = +00 and £ = e 11 to x = —00. 
The complex potential for the non-symmetric case is 



W(C) = *7/2 + log 



where 

e(C) = (i-C) n (i-? 2m C)(i-<? 2m /0 = 



+ 00 



e(e-^c)e(aC) 

e(ge-*rO 9(C) ' 

tf 4 (iog(VCM9) 



and d 4 (w,q)= ^] (-l)"g n2 exp(2 



nw 



n— — oo 



is the Jacobi theta function [16( ■ It is easy to verify that 
W(C) satisfies the boundary conditions indicated above. 
Finite-parametric solutions in this case have the form 



z{t,4>) = r(f) + »7(*)/2 + log 



e(e*(*-7(t))) 



0( e »0) 

Y^k lo g 6(a fc (i)e-^) + a k lo g e(a fc (f)e*)]. (10) 



fc=i 



Here all |ofc| < 1, and J2k=i a k = 0- Since y is multiple 
of 7r, when |£| = 1, then r is purely real. Integrating the 
equations of motion, we obtain N + 2 complex constants 
of motion: j3 k = f(t, q/a k ) for k = 1, N; P + = f(t, q), 
and f3- = /(i,ge* 7 ), where 



Pk 



i 1 /2 + \og[Q{qe^a k )/Q{qa k )}- 



N 



^ [a; log Q(aia k /q) + m log 9(ga;/a fc )] , (11a) 



i=i 



p± = t - 2t + i 1 /2 ± log [6(g e*)/e(g)] 



jv 



^ [aj lo g e(e- l T± ai/q) + a, log 9(ge^±a ( )] (lib) 



1=1 
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FIG. 3. Examples of bubble evolution: (a) symmetric solution 
and (b) asymmetric solution. 

where 7+ = 1 and 7_ = 7. The constants, j3+ and /?_ 
are not independent, since Zm/3+ = Im/3_ = 7/2 + 
Im 2/=i a ' 1°§ a (- The formulae (11) constitute the full 
dynamics of a^, 7, q, and r. The bubble area, A, while 
fixed in time, is not an independent constant of motion; 
it is neatly expressed through other constants as 

jv 

A/ ir = p+-p- + 2TZe ^kh- 
fc=i 

The attractor. After eliminating r, by subtracting (3± 
from we see at once from (11) that all ^ q move 
to the point qe n , when t — > 00. The pftoints q and 
qe n are the repeller and the attractor frespectively for 
the dynamical system expressed by (11). Taken at the 
limit, t — > 00, these two points are the only fixed points 
of this system. If one of o^'s, say ai, was initially at 
the repeller, a\ = q, then /3i diverges, and U = 2/(1 + 
Qi) 7^ 2, when t — > 00. After relocating a\ out of q, the 
asymptotic velocity of the bubble eventually reaches the 
selected value, U = 2. 

The solution (10) describes symmetric bubbles, when 



/3 = 7T, thus (10) contains (6), albeit in a different formu- 
lation. A symmetric bubble evolution is shown in Fig.^, 
where the asymptotic shape corresponds to the Taylor- 
Saffman bubble with U = 2, as described by (8). In 
Fig. [3Jd we show an asymmetric solution whose steady 
shape (reached in the long-time limit) corresponds to the 
Tanveer asymmetric bubble [||,[l2j with U — 2. 
Final remarks: 

(i) Since all move toward the same point, ak — > qe 11 , 
for t — > 00, the sum in (11) vanishes, because X)fc=i a k = 
0, so j(t) —> Imf3 + = Im[3-, which, being a constant of 
motion, can be any number chosen freely at t = 0. Thus, 
it remains to be seen how to centralize a bubble in our 
framework (if possible at all), so that 7 — > ir as t — > 00. 

(ii) We do not present here the U = 2 selection for 
the asymmetric bubble with infinitely-parametric non- 
singular solutions, that is the analog of (9) but with use of 
theta-functions. While the algebra is lengthy, the result 
is the same, U = 2, in the limit t — > 00, as expected. 

(iii) Presented results (and the result of [U , which is a 
singular limit of this work: v — > 1 + e, for < f « 1) 
unambiguously indicate that stability of the selected pat- 
tern, with respect to the rest of the family, is built in the 
Laplacian growth without surface tension. Surface ten- 
sion is just one of infinitely many perturbations (perhaps 
the most relevant), which kicks the system toward the 
attractor, while also regularizing high curvatures. 

(iv) The selected pattern, while linearly unstable, is 
stable asymptotically: if perturbed, it eventually recovers 
its original shape. It is tempting to regularize our theory 
(by means of surface tension or otherwise) and to see 
which subset of the presented solutions will survive. 

One of the authors (M.M-W) would like to thank the 
Max Planck Institute for Complex Systems (Dresden) for 
hospitality during preparation of the part of this work. 
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